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Abstract 

We review the recent advances in the Hamiltonian formnlation of lattice 
gauge theory for approaching the continuum physics. In particular, vacuum 
wave function and glueball spectrum calculations by coupled cluster method 
with truncation scheme preserving the continuum behavior are described. 


1 Introduction 

Lattice gauge theory (LGT) has now developed into one of the most informa¬ 
tive method for non-perturbative aspects of strong interactions. It is a funda¬ 
mental theory based on hrst principle of strong interactions. In 1974, Wilson 
proposed the Lagrangian formulation of LGT where the gauge held the¬ 
ory is discretized on the D-dimensional space-time lattice. In 1975, Kogut and 
Susskind derived a lattice Hamiltonian in which only the (D-l)-dimensional 
space is discretized while the time variable remains continuous. Theoretically, 
these two approaches are equivalent in the continuum limit. 

For simplicity, we consider only the case of pure gauge helds. In the con¬ 
tinuum, the action is 


s=-\j d-^x ( 1 ) 

Gontinuum quantum held theory with such an action suhers from divergence 
problem, which has to be regularized and renormalized. In perturbative QGD, 
the theory is regularized at high momenta, and renormalization group equation 
gives the ehective coupling constant ^ as a function of the momentum scale 
/i, i.e. g = gi^p). For SU(N) group, g goes to zero when /i —> cx), so that the 
high energy strong interactions can be well described by perturbative QGD. 
For low energy phenomena, such as quark conhnement and glueball spectrum, 
where g becomes large, non-perturbative methods must be employed. 
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LGT uses the space-time discretization to regularize the continuum theory. 
On the lattice, the gauge held A is replaced by the link variable Ui{xo) = 
exp{ig dxiAi), where a is the lattice spacing. Renormalization group 
equation predicts the so-called scaling relation g = g{a). g goes to 0 as a —0. 
Therefore, the continuum limit corresponds to the weak coupling limit of the 
lattice theory. In this limit, the evolution of g{a) is the same as the evolution 
of g{fi) in the continuum theory, with a ~ /i“^. For a 4D SU(N) theory, the 
scaling relation is 


_1 

Aia = g ^‘'oe {1 + cig'^ + ...). (2) 

The constants Bq, bi and Cj can be calculated by continuum perturbation theory. 
Al is the lattice scale parameter with dimension of mass, and it can be related 
to the continuum scale parameter such as When a varies, the lattice 

coupling constant g must vary accordingly in order to give the same physical 
results for continuum physics. On the lattice, one calculates dimensionless 
quantities such as am, with m some physical mass. In the weak coupling 
region, am should scale as const, x A^a^g). Correct scaling behaviors is crucial 
for extracting physical results from lattice calculations. 

The last two decades has seen a lot of progress in Monte Carlo (MC) 
simulations of LCT in the Lagrangian formulation, while analytic and numer¬ 
ical investigations of Hamiltonian formulation were not so active. Nowadays, 
MC simulation of lattice QCD is able to reach (3 = Q/g'^ k, 6.4, but it is 
unfortunately not weak enough for the scaling behavior to be satished, and 
improvement is a pressing task. In recent years there have been some remark¬ 
able advances in the Hamiltonian formulation for approaching the continuum 
physics. It now seems that the Hamiltonian formulation deserves much atten¬ 
tion. 


2 Hamiltonian versus Lagrangian Formulation 

The Lagrangian formulation starts from the path integral Z = 

Wilson [|I| proposed an action 

S=^T.Tr(Uv + Ul-'^)^ ( 3 ) 

y p 
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where Up is the ordered product of U around an elementary plaquette. The 
mass spectrum can be obtained by computing the correlation function C{t) = 
(<h(t)<h(0)), where <h(t) is some appropriate operator with nonzero projection 
onto the lowest excited state. Then 


C{t) = Z-^ = Tr[e^'<h(0)e-^*<h(0)] 



( 4 ) 


n 


For to hold at not too large t, <F must have sufficiently large projection 
onto the lowest excited state. Furthermore, statistics must be sufficient and 
hnite size effects must be under control. Despite existing problems, the La- 
grangian formulation simulated by MC algorithm is still the most efficient way 
for obtaining low energy hadron mass spectrum. 

In contrast, the Hamiltonian formulation starts from the lattice Hamilto¬ 
nian [H 



( 5 ) 


where Ef is chromo-electric helds on the I link, and it is also generator of the 
gauge group. The mass spectrum can obtained directly by solving the eigen- 
equation H'^[U] = e^T[[/]. Here is the eigenvalue of H. When a —> 0, a 
huge number of gauge conhgurations are correlated, and it is very difficult to 
diagonalize the Hamiltonian with sufficient accuracy. For this reason, there 
is no satisfactory (3-1-1 )-D result in this formulation up to now. In a feasible 
calculation, the gauge conhguration space has to be truncated. An inappro¬ 
priate truncation scheme often violates the continuum limit, and destroys the 
scaling behavior for the physical quantities. Therefore, special care must be 
taken when choosing a truncation scheme. 

On the other hand, the Hamiltonian method does have some advantages 
over the Lagrangian method. It is relatively simple to obtain wave functions 
of the hadronic states, and it is likely that more physical information on the 
glueballs can be derived in this formulation. 
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3 Brief Comments on Earlier Investigations 

We would like to mention here some earlier work on the Hamiltonian formu¬ 
lation: 

a) Strong coupling expansion with Fade approximants 0. Since the strong 
coupling expansion diverges at small g, some approximants must be used to 
extrapolate the strong coupling series to the weak coupling region. This intro¬ 
duces uncertainties and no conclusive results have been obtained. 

b) Eigen-equation method with truncated gauge conhgurations proposed by 
Greensite [Q. This is just the coupled cluster method with a special truncation 
scheme. In |^, only the strong coupling region was explored. However, the 
scaling region was not reached and no concrete results were obtained. 

c) Variational method. The vacuum energy eo is obtained by minimizing the 
expectation value of H in the trial vacuum state |f2), and the mass gap is 
obtained by minimizing the expectation value of H in the trial excited state 
I'll) with (HjT) = 0. Some earlier results [Q were consistent with asymptotic 
scaling predictions. However, since the variation energies depend strongly on 
the choice of trial wave functions, it is not clear whether this method can give 
reliable results for the mass spectrum in a systematic way. 

d) Models with exact ground state . Some modihed Hamiltonians with clas¬ 

sical continuum limit the same as that of the Kogut-Susskind (K-S) Hamil¬ 
tonian were proposed. These models possesses an exact ground state of the 
form |f2) = where |0) is the fluxless state and S{U) is some gauge 

invariant operator. The mass gaps were obtained by minimizing the excitation 
energies, which was shown to have good scaling behavior in (2-|-l)-D theories. 
The problem is that the modihed Hamiltonian differs at quantum level from 
the K-S Hamiltonian by a relevant operator and hence belongs to a universal 
class different from the K-S Hamiltonian. 

e) Linked cluster expansion method [0]. While the Hamiltonian can not be 

diagonalized exactly, its sub-matrix on a hnite set of strong coupling 

basis Tj can be diagonalized. This is a variant of strong coupling expansion, 
and some approximants must be used to extrapolate the results into the weak 
coupling region. However, such an extrapolation might lead to uncertainties. 

f) Unitary transformation and variational method for LGT with fermions [Q. 

A variational form of hadronic wave functions that takes into account the 
effect of sea quarks was proposed |tk) = Meson spectrum and chiral 

condensates were obtained in low dimensional cases. They were consistent with 
scaling predictions. This is a hrst step towards the more systematic coupled 
cluster method. 
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To summarize, much efforts have been made in earlier investigations, but it 
seems difficult to use them to study the continuum physics of a realistic theory 
such as QCD 4 . More systematic and unambiguous methods are required. 

4 Recent Advances: Coupled-Cluster Method 

In the 90’s, one has seen considerable developments in Hamiltonian LGT for 
approaching the scaling region. One of the most promising analytical method 
is the coupled-cluster method with an appropriate truncation scheme [^. We 
will discuss this method is some details. 

4.1 Coupled-cluster method (CCM) 

Let the vacuum state be 


|Ji) = e«PI|0), (6) 

The eigen-equation H\Q) = en|0) becomes 

2 2(7 

i:{ii5r. [£,".«!] + [.Er, filler. fli} - -ri:fr(E,+ = -j£„. ( 7 ) 

i 9 p 9 

R[U] is composed of various gauge invariant Wilson loops, and it can be ex¬ 
panded in a series of graphs 


11 , 12 


R[U] = ^ ^ a, ^ GnAU]. (8) 

n n,i X 

For example, in (2-|-l)D SU(2) theory, Gi^i = □, ^ 2,1 = ^ 2,2 = | | ,etc. 

The eigen-equation (|^ is a system of nonlinear equations for the coefficients 
Gn,i- In practice, this equation must be truncated at some finite n. Inappropri¬ 
ate truncation scheme would violate the continuum limit of (13) and destroy 
the scaling behavior. 

4.2 The continuum limit of a graph 

The continuum limit of a generic graph Gn,i has the form 
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G„,i\U] = gV[An,i Tr{r-) + Tr{Vrf + ...]. 


(9) 


For example, the elementary plaquette behaves as 


Up = ReTr{PU^o = ReTr[l + i f gA ■ dx - hf gA ■ dxf + ...]. (10) 

Ju 2 Ju 

Note that up to the A^ term, it is not necessary for path to be ordering because 
of the trace. Once the A^ term is determined, we can supply the A^ and A^ 
terms to make up gauge invariant expressions. Thus with xq at its center. 


A ■ dx 



dxpdxi, [dpA^{xQ + x) 


d^Ap{xQ + x)] 


= + VDRp^xo) + 0(a^ 


( 11 ) 


Therefore, (pTf) becomes 


Tr{U, - 1] = + ... (12) 

The continuum limit of other graphs can be calculated in a similar way. There¬ 
fore, the operator R has the continuum limit 


R[U] = -jd°-^x + ...]. 


(13) 


The vacuum state for the SU(2) gauge theory was hrst investigated by 
Greensite and later by Arisue [|^ using simulation method. Nice scaling 
behavior for po and /i 2 was obtained in the 3D theory |^. These coefficients 
approach constant values in the scaling region. 


4.3 Truncation scheme preserving the continuum limit 

The eigen-equation (J^) has to be solved in some truncation scheme. If in the 
expansion (|^), we let n be the order of a graph , then we have 
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( 14 ) 


[E^, Rn]] G Rn + lower orders terms, 

i 

, R\[E°, R] G Rni+n 2 + lower order terms. (15) 

i 

Thus the last term must be truncated such that the orders of graphs appearing 
in (^ do not exceed some hnite order N. A conventional prescription [^, p!3[| 
was that in ([TED , only graphs with order N were preserved, while graphs with 
order greater than N were discarded. 

We El have shown that these truncation schemes violate the continuum 
limit of the term (pAj) , and hence destroy the scaling behavior of po and p 2 • The 
essential point is that, since a general graph Gi has the continuum limit (§), 
the continuum limit of [Ef ,Gi\[Ef',Gj\ is 


GJ [£■,“. Gj] K Tr{V}=^_,Y + ... 


(16) 


For this equation to be valid, all graphs generated from the l.h.s. must be 
included. On the contrary, if some graphs were kept and others were discarded, 
the continuum limit would change into a^g‘^Tr{E^yE^y), leading to seriously 
wrong continuum behavior. 


Based on this consideration, we proposed a truncation scheme that 
respects the continuum limit. 


N 


ll{\Er^Er,Y^R„(U)]]+ Y. |G?,fl„.(G)l[Br.-R™(G)]} 

I n ni-\-n2<N 


- (1'^) 
y p y 

To conhrm the validity of this method, we must test the convergency of the 
truncation scheme at large N and check the scaling behavior of various quan¬ 
tities, in particular, the coefficients po, and /i 2 of the vacuum wave function 
and the glueball masses m(J^^). 
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4.4 Results 


We have examined this method for several lattice field systems: 

a) ( 2 +l)D U(l) theory [|T^]. hq and the mass gap niA are calculated up to 8 th 
order. The results show clear tendency of convergence and exponential scaling. 
The mass gap agrees with MC and other analytic results. (See Fig. 1). 

b) (2+l)D SU(2) theory [|^, |^. Reasonable results with power scaling for 
/io, /i2 and ms are obtained even at the 3rd order. Convergent results are 
obtained up to 7th order. The continuum limit of ms ~ 1.4 is somewhat 
lower than Teper’s MC results. (See Fig. 2). 

c) (2+l)D SU(3) theory [T^. Up to the 3rd order, we obtained j^o, m(0’''+) 

and m(0 ), consistent with power scaling, m(0’''^) ~ 2.09e^, and m(0 ) ~ 

3.71e^. (See Figs. 3, 4). 

d) (1+1)D 0(3) a model |]^. Up to the 8 th order, the triplet mass mx 

is consistent with MC results, indicating that asymptotic scaling occurs at 
rather large values of At smaller values of /?, there exists a scaling region 
according to the whole (3 function. We conhrm that oc const, in this 

region. (See Fig. 5). 

e) (1+1)D 0(2) model fl^ (Hamiltonian formulation of 2D XY model). Con- 
vergency was examined up to the 10th order. There is a clear signal of KT 
phase transition at 5 ^^ ~ 1. The critical point and critical exponent are consis¬ 
tent with results by other methods. (See Fig. 6 ). 

f) Estimates of QCD 4 glueball masses from QCD 3 results by dimensional re¬ 
duction |l^. In the strong coupling limit or large W limit, the hxed time 


vacuum expectation value of an operator 0{U) in D dimensional conhnement 
theory (2 < D < 4) correponds to the path integral expression < 0{U) > 
in (D-l)-dimensional theory. For long wavelength excitations, the mass ratios 
should be approximately the same for (3-|-l)D and (2-|-l)D. Combining our 
updated QCD 3 data m(0’''+) ~ (2.15 ± 0.06)e^, and the recent MC data for 
the QCD 3 string tension a = (0.554 ± 0.004)e^, we obtain 


V ^ ^ 3.88 ±0.11, 18 

Va 

which is to be compared with the recent MC data for QCD 4 : m{0 ~^~^)/~ 
3.95. We see that the dimensional reduction relation works quite well in 
this case. If we further use the MC result for the QCD4 string tension a ~ 
0.44 GeU, we obtain for QCD 4 : 


m(0++) 1.71 ±0.05 Gev, 


(19) 












which is consistent with the IBM data = 1.740 ± 0.071 [pT| . 

g) Preliminary calcnlations of QCD spectrnm in (3+l)D [^. We have also 
calcnlated the 0 ’''+, 0 and l+“ glneball masses in QCD 4 in low order ap¬ 
proximations. Althongh the scaling behavior for these glneball masses has not 
been achieved, we do observe plateans for glneball mass ratios in /3 G (6.0, 6.4). 
From the plateans we obtain 

M(0—) M(l+-) , , 

’ = 2.44 ±0.05 ±0.20, —= 1.91 ± 0.05 ± 0.12, 20 

M(0++) ’ M(0++) ’ ^ ^ 

where the hrst error is the error of the data in the platean, the second error 
is the estimated error dne to hnite order trnncation. Onr resnlts are in good 
agreement with the MC and t-expansion resnlts. 


5 Other Developments 


a) Improved Hamiltonians. The idea of improved lattice actions was proposed 
in the early 80’s 1^^. Recently the improved actions have attracted mnch 
attention [^. The pnrpose is to pnsh the 0(a) errors to higher orders and 


remove the tadpole lattice artifacts. Using improved actions, one can obtain 
resnlts on coarser lattice, which was previonsly obtained on hner lattice. The 
compntation time may then be greatly rednced. Starting from Lepage’s im¬ 
proved Lagrangian, an improved Hamiltonian [Q can be straightforwardly 
derived nsing the transfer matrix method or the Legendre transform method. 
However, the color electric energy becomes an inhnite series with long range 
terms. This dehciency can be cnred by the Legendre transformation of a snit- 
able improved Lagrangian with inhnite time-like terms, yielding the same order 
of improvement bnt with only local and nearest neighbor interactions 
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b) Universality of LGT in Hamiltonian and Lagrangian formnlations p5 
Hamer et al. recently derived the relation between the conpling constants 
qh and ql in both formnlations. Relations between velocity of light in both 
formnlations were also derived. Using these relations, a comparison was made 
for (2-|-l)D SU(2) theory between varions Hamiltonian calcnlations and MC 
resnlts in the Lagrangian formnlation. A striking demonstration of nniversality 
between both formulation was obtained. 

c) Calculation of hadronic structure functions 


ZB|. A major difficulty 


m 


Hamiltonian LCT is that a great number of correlated conhgurations are in¬ 
volved. Recently Kroger and Scheu showed in a scalar model that by using 
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specific reference frame (e.g. Breit frame), the number of relevant configu¬ 
rations may be greatly reduced, so that they could obtain nice results for 
distribution function consistent with the theoretical prediction. 


6 Conclusions and Outlook 

Recent work on Hamiltonian LGT shows that with appropriate calculation 
schemes (CCM with continuum limit preserved truncation, Breit frame, etc.), 
we can efficiently enter the scaling region and obtain reliable results. More 
calculations on (3-|-l)D should be performed for further establishing the ad¬ 
vantage of the Hamiltonian formulation in spectrum and wave function cal¬ 
culations. As the number of correlated, conhgurations increase considerably 
from (2-|-l)D to (3-|-l)D, improved Hamiltonians may be necessary for efficient 
calculations in (3-|-l)D. 
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This work was support by National Natural Science Foundation under grant 
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mittee under grant number jiao-wai-si-liu [1996] 644. 


References 

[1] K. Wilson, Phys. Rev. DIO (1974) 2445. 

[2] J. Kogut, L. Susskind, Phys. Rev. Dll (1975) 395. 

[3] T. Banks et al., Phys. Rev. D15 (1977) 1111. 

[4] J. Creensite, Nucl. Phys. B166 (1980) 113. 

[5] H. Arisue et ah. Prog. Theor. Phys. 70 (1983) 229. 

[6] S.Cuo et ah, F%s.i?et'.D38(1988)2571;F/i?/s. Rev. D44(1991)1296. 

[7] C. Hamer, A. Irving, T. Preece, Nucl. Phys. B270(1986) 553. 

[8] Q. Chen, X.Q. Luo, Phys. Rev. D42 (1990) 1293; X.Q. Luo et ah, Phys. 
Rev. D50 (1994) 501. 

[9] J. Creensite, Phys.Lett. B191(1987)431. 


10 



[10] H. Arisue, Phys. Lett. B280 (1992) 85. 

[11] S. Guo et al., Phys.Rev. D49(1994)507; Nucl.Phys.'B(P.S.)47{1996)827. 

[12] C. Hamer, W. Zheng and D. Schiitte, |hep-lat/9603026| . 

[13] C. Llewellyn-Smith, N. Watson, Phys.Lett. B302(1993)463; R. Bishop et 
ah, Phys.Rev. 048(1993)887. 

[14] X. Fang, J. Liu, S.H.Guo, Phys.Rev. 053(1996)1523. 

[15] Q.Ghen et al.,P%s.i?ev.050(1994)3564; Fang, Guo, Luo, in progress. 

[16] Q. Ghen, X.Q.Luo, S.H.Guo, Phys.Lett. B341(1995)349; Phys.Lett. B348 
(1995) 560; Nucl.Phys.B(P.S.)47 (1996) 274. 

[17] X. Fang, S.H.Guo, J.M.Liu, Phys.Rev.T)b4 (1996) 6521. 

[18] J.M.Liu, X.Y.Fang, S.H.Guo, Phys. Rev. B (1997) 

[19] X.Q.Luo, Q.Z.Ghen, Mod. Phys. Lett. All (1996) 2435. 

[20] X.Q.Luo et ah, Nucl. Phys. B(P.S.)53 (1997) 243. 

[21] J.Sexton, A.Vaccarino, D.Weingaren, Phys. Rev. Lett. 75 (1995) 4563. 

[22] K.Symmanzik, Nucl. Phys. B226(1983)187, 205. H. Hamber, G.M.Wu, 
Phys.Lett. B133 (1983) 251. 

[23] G.P.Lepage, 1996 Schladming Winter School Lectures. 

[24] X.Q.Luo, S. Guo, H. Kroger, D. Schiitte, Nucl. Phys. B(P. S.), to appear. 

[25] G. Hamer et ah, Phi/s.Ren.054(1996)2395; Phys.Rev.D53{1996)7316. 

[26] H. Kroger and N. Scheu, Nucl.Phys.'B{P.S.)53{1997)321. 


Figure Captions 

Fig.l. Mass gap tha in (2+l)D U(l) theory. 

Fig.2. 0~^ glueball mass in (2+l)D SU(2) theory. 

Fig.3. Goefficients in QGD 3 vacuum wave function. A; strong coupling 
expansion; x: third order; Dot lines: mean values in the scaling region. 
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Fig.4. QCD3 glueball masses marked by x. Dot lines: mean values in the 
scaling region; Dash lines: strong coupling expansion; A: third order from 
a set of connected graphs. 

Fig.5. Triplet mass gap rriT in (1+1)D 0(3) nonlinear model. 

Fig.6. Mass gap uia in (1+1)D 0(2) model. Solid circles: fit to the KT 
form. 
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This figure "fig 1-1.png" is available in "png" format from: 
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